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FOREWORD 


This  interim  technical  report  was  prepared  by  Electro-Optical 
Systems,  Inc.,  Pasadena,  California,  on  Contract  AF  33(657)-7940  for 
the  Aerospace  Research  Laboratories,  Office  of  Aerospace  Research, 
United  States  Air  Force.  The  research  reported  herein  was  accomplished 
on  Task  7063-03,  "Energy  Exchange  Phenomena  in  Electric  Arc  Discharges" 
of  Project  7063,  "Mechanics  of  Flight"  under  the  technical  cognizance 
of  Mr.  Charles  A.  Davies  of  the  Thermo-Mechanics  Research  Laboratory 
of  ARL. 


ABSTRACT 


Simple  analytic  expressions  are  derived  for  describing  the 

radiation  from  an  arc  confined  in  a  cylindrical  constrictor  and  operating 

in  the  general  "gas-heating"  or  "inlet"  regime.  These  expressions  were 

obtained  by  postulating  linear  relationships  among  the  plasma  properties 

a,  u  and  (|),  the  electrical  conductivity,  the  power  radiated  per  unit 

T 

volume,  and  the  thermal  conductivity  integral  ^  Kdt,  respectively. 

o 

For  initial  comparisons,  radiation  from  "fully  developed"  arcs  in 
argon  and  nitrogen  has  been  measured  as  a  function  of  the  arc  current 
and  the  gas  pressure  in  the  confining  cylinder.  These  measurements  are 
compared  to  the  appropriate  limit  in  the  analytical  predictions  and  the 
correlation  is  found  to  be  good,  except  for  the  electric  field  dependence 
upon  current.  It  has  been  found  possible  to  improve  the  correlation  by 
modifying  the  results  of  the  analysis  to  include,  in  an  approximate 
manner,  the  nonlinear  relationships  that  actually  exist  among  the 
properties  of  the  plasma.  Similarity  parameters,  which  are  inferred 
from  the  analysis  and  approximately  derived  from  the  basic  equations, 
are  applied  successfully  to  the  measurements  with  argon  and  nitrogen 
and  to  extensive  calculations  with  nitrogen.  Radiation  measurements  for 
comparison  with  the  complete  analysis  remain  to  be  accomplished  in  the 
"gas  heating"  or  "inlet"  regime,  but  the  work  to  date  offers  hope  that 
the  analysis  has  included  the  main  features  of  the  arc-gas  interaction. 
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INTRODUCTION 


The  new  technologies  of  hypersonic  flight,  space  flight,  rocketry 
and  light  sources,  as  well  as  some  aspects  of  nuclear  fission  and  fusion 
work,  have  created  several  new  requirements  for  sources  of  hot  gas  or 
plasma  above  the  combustion  temperature  regime.  Thus,  the  heating  of 
various  gases  by  means  of  electric  arcs,  typically  to  temperatures 
between  3000  °K  and  50,000  or  higher,  at  pressures  of  from  1  atmos¬ 
phere  to  several  hundred  atmospheres,  h.as  become  an  important  technology, 
resulting  in  many  specialized  forms  of  arc  plasma  sources,  or  arejet 
devices.  At  these  pressures  and  temperatures,  the  power  that  is  radi¬ 
ated  per  unit  volume  from  t!ie  gas  can  become  the  major  component  of  heat 
transfer  in  the  jas.  Techniques  must  be  devised  for  properly  incorpor¬ 
ating  this  phenomenon  into  the  design  concepts  of  arc  heater  and  arc 
light  source  devices. 

Considerable  theoretical  work  lias  been  done  in  an  attempt  to  com¬ 
pute  the  power  radiated  per  unit  volume  from  tlie  gas  as  a  function  of 

(1  2) 

the  gas  temperature  and  pressure'  ’  .  There  is  general  agreement  that 

the  major  portion  of  the  radiation  results  from  a  combination  of  two 
types  of  radiative  transitions  of  "free"  plasma  electrons.  The  well- 
known  bremsst rah  lung  continuum  results  from  electrons  being  deflected 
by  coulomb  interactions  with  positive  ions  and  transforming  kinetic 
energy  into  radiation;  added  to  this  radiation  is  that  due  to  radiative 

Manuscript  released  by  the  authors  3  June  196A  for  publication  as  an  ARh 
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recombinations  of  electrons  into  capture  orbits.  An  equation  describing 
these  two  processes  and  which  allows  computation  of  the  power  radiated 
by  the  gas  in  terms  of  known  physical  constants  has  been  derived  by 
Unsold^^^  using  the  classical  Kramers  theory.  Yos^^^  has  furtlicr  cor¬ 
rected  this  formula  to  match  very  limited  experimental  data,  and  has 
computed  and  plotted  the  power  radiated  per  unit  volume  for  a  number  of 
gases  as  a  function  of  temperature  and  pressure. 

Because  of  the  rather  approximate  nature  of  the  ca Iculations,  which 
may  well  allow  an  error  in  the  absolute  magnitude  of  a  factor  of  3  or 
greater,  it  is  highly  desirable  to  carry  out  experiments  involving  the 
measurement  of  radiation.  It  haS  been  found  that,  in  practice,  experi¬ 
ments  are  very  difficult,  and  progress  has  been  quite  slow  in  finding 
techniques  that  will  allow  a  simple  one-to-o-  comparison  between  theory 
and  experiment  including  radiation.  In  this  report,  one  approach,  with 
some  preliminary  results,  is  discussed  in  detail.  Here  the  radiation  is 
from  argon  and  nitrogen,  heated  by  an  arc  in  a  confining  cylinder.  The 
pressure  and  cylinder  diameter  are  maintained  at  values  such  that  the 
gas  should  be  "optically  thin,"  i.e.,  a  negligible  amount  of  self¬ 
absorption  occurs.  Analytical  and  experimental  results  are  then 
compared  on  the  basis  of  parameters  derived  from  linearized  approximations 
in  the  theory  and  from  similarity  considerations. 
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APPROACH 


The  configuration  that  has  been  chosen  for  the  detailed  study  of 
the  radiation  from  the  gas  is  shown  schematically  in  Fig,  1*.  In  the 
correlation  between  theory  and  experiment^  it  is  assumed  that  the  gas 
flow  rate  through  the  cylinder  is  low  enough  so  that  all  of  the  power 
dissipated  in  the  discharge  is  transferred  to  the  segment  walls  by 
conduction  and  radiation.  However^  in  formulating  the  analysis  it  is 
desirable  to  be  able  also  to  describe  the  region  where  the  gas  is 
being  heated^  so  that  eventually  the  parameters  of  the  heater,  e,  g.,, 
constrictor  length  and  diameter,  can  be  specified.  This  requires  that 
an  attempt  be  made  to  solve  the  steady-state  energy  equation  with  vari¬ 
ation  in  both  the  axial  and  radial  dimensions.  In  order  to  obtain 
simple  analytical  solutions  which  are  very  desirable,  linear  relations 
are  postulated  among  the  quantities  I,  C,  and  u,  where  (|>  is  the  thermal 
conductivity  integral,  c  is  the  electrical  conductivity,  and  u  is  the 
power  radiated  per  unit  volume  of  gas.  Such  linear  relations  can  be 
expected  to  hold  only  over  very  limited  ranges  of  temperature. 

However,  useful  relations  and  parameters  may  be  obtained  from  the  line¬ 
arized  solution  which  migh.t  then  suggest  the  formulation  of  more  general 


This  classical  "Maecker-type"^^^  column,  and  the  measurements  and 
theory  associated  with  it  are  detailed  or  referenced.'*''  ' 
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similarity  parameters  that  have  a  much  wider  range  of  application  than 
would  the  linear  solution  itself.  Another  point  uf  view  can  also  be 
adopted  such  that  the  linear  relations  are  only  average  values  of  the 
essentially  nonlinear  relations  among  the  variables.  Once  the  solution 
is  obtained^  these  average  values  can  then  be  replaced  by  a  more  accurat 
nonlinear  relation  in  the  solution.  The  validity  of  this  procedure  can 
best  be  checked  by  self-consistent  comparisons  between  the  experimental 
results  and  the  "modified"  solution.  Tliis  comparison  is  carried  out  in 
the  following  sections. 
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LINEARIZED  ANALYSIS 


An  approximr'-c  solution  is  known  for  tne  axial  and  radial  distri¬ 
bution  of  gas  properties  in  the  steady  flow  of  a  gas  in  a  cylinder  with 
an  axial  electric  arc.^^^  The  theory  will  be  expanded  to  include  radia¬ 
tion  in  the  gas. 

The  cntl.alpy  equation  states  tliat^  fui  steady  fiow^  the  net  rate 
of  heat  outflow  from  a  unit  volume  by  conduction,  convection  and 
radiation  is  equal  to  tlie  power  input  due  to  ohmic  heating.  If  radial 
convection  and  axial  conduction  are  neglected,  the  equation  in  cylindrical 
coordinates  is: 


■57  + 


+  u 


.2 


(n 


where 


T 


^  =  thermal  conductivity  integral,  “  j 

o 

u  =  power  radiated  per  unit  volume  (self-absorption  neglected) 
The  following  approximate  solution  is  valid  for  the  "gas  heating" 
or  "entrance"  region  for  which  the  arc  (conducting  region)  does  not  fill 
the  cylinder. 

Assume  that  5w  is  constant,  which  implies  the  neglecting  of  radial 


flow,  Assum.e  that 


where  is  a  function  of  pressure  only. 

\  y  /d\  f  ^ 

Assume  that  c  is  zero  for  i  less  than  A,  and  a  -  (tt)  I  '  fA  - 

h\  WAPa/ 

for  greater  than  where  y^)  i®  ^  constant  evaluated  at  atmospheric 
pressure  (for  convenience). 

By  using  the  fact  that  the  integral  of  the  current  density  j  over 
the  cross  section  is  the  total  current  I,  and  assuming  that  the  electric 
field  E  =  j/a  is  independent  of  radial  position^  we  can  write: 


Al’  of  the  above  assumptions  are  essentially  the  same  as  those  made  by 
Stine  and  Watson  except  for  the  explicit  inclusion  of  the  pressure 
dependence  of  c\  the  pressure  dependence  is  likely  to  be  small  but  is 
included  for  later  reference. 

The  radiation  u  in  Eq,  (1)  is  the  net  power  radiated  per  unit  volume, 
Assume  that  the  column  is  optically  thin^  and  hence  that  the  self¬ 
absorption  by  the  gas  can  be  neglected.  We  assume  that  the  total 
specific  radiation^  u,  is  zero  for  less  than  and 


“  ‘wJk 

for  ^  greater  than  where  Is  a  constant  evaluated  (again^  for 

convenience  only)  at  atmospheric  pressure.  The  value  of  in  the 

above  assumptions  is  obtained  by  extrapolating  the  straight  line 
approximations  between  c  and  and  between  u  and  (at  atmospheric  pressure 
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to  the  I  axis  and  adjusting  them  so  that  the  same  value  of  is  obtained 

from  the  two  lines.  The  value  of  r(-  r,)  corresponding  to  then 

defines  an  "effective  arc  radius."  Actual  values  for  c('|>)  and  u(;|i)  may 
be  obtained,  for  example,  from  the  data  given  or  quoted^^^  for 

various  gases. 


With  the  above  assumptions,  Eq.  1  becomes: 


for  r  s  and 


1 

r  Sr 


(3) 


for  r  >  rj^. 


The  solution  of  Eq.  (2)  for  ({>  =  (ji^^  at  r  =  r^^  and  z  =  0  is 


I  (2.4  r/r^)  f  (z) 


(^0 


where 


Notice  that  even  when  radiation  is  included,  the  solution  takes  the  form 
of  a  product  solution;  tlius  the  solution  at  any  radius  is  tlie  product  of 
a  function  of  r  witli  a  function  of  z. 
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Eq,  (5)  may  be  written  in  the  following  form: 

c/  N  /I  11-5  z/z 

f(z)  >=  (1  -  e  '  o) 


where 


If  we  write  (pw)  =  - -  where  m  is  the  mass  flow  rate^ 

n 


(6) 


(7) 


then 


(8) 


In  brief,  the  axia  1  depei.dence  of  the  various  quantities  will  remain 
the  same  as  in  the  original  calculations  of  Stine-Watson^^^  provided 
that  the  quantity  is  now  g’ven  by  Eqs.  (7),  (8).  As  a  result,  it  will 
suffice  to  consider  only  the  radial  dependence  in  the  following  discussion. 

The  complete  radial  dependence  of  (i.e,,  the  asymptotic  value 
for  large  z)  is 
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(Hi) 


1  J  (2. A  r/r,) 


2n  J^(2.4) 


\2.u)  vv;,.,pj 


1/2 


r  &  r. 


(9) 


or 


■2,4  1  In  (r/tj^) 


2T7r 


1  + 


I/; 


<  r  i  R 


(10) 


where  the  "outer  solution"  (no  ohmic  heating  or  radiation)  was  found  by 
setting  d^/ 5z  ==  0  and  matching  &|i/ ?r  with  the  "inner  solution." 

To  initially  compare  the  predictions  of  the  linear  theory  with 
experimental  measurements^  only  the  asymptotic  arc  column  region  with 
the  solution  given  in  Eqs.  (9j  10)  will  be  considered.  (No  experimental 
data  including  radiation  measurements  are  available  for  the  inlet  region 
as  yet.) 

For  later  reference,  the  complete  predictions  of  the  theory  applied 
to  the  asymptotic  arc  column  regime  are  listed  below: 

Ratio  of  total  power  radiated  per  unit  lengtii  to  total  power  conduc t ed 
per  unit  length,  § 


5  = 


(H) 


Total  power  radiated  per  unit  length,  P 

^ rfl  (j 


(12) 
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J^(2.4  r/r^) 


2rj^(2.4) 


( 


10 


"Radius  of  arc," 


R/r^  =  exp 


fjj  J 


(19) 


=  exp 


[ 


2n  (fj  (1  +  5) 


El 


] 


(20) 


Average  gas  enthalpy,  h 


h  -h, 
avg  1 


(t)  2.4 


TI 


/daW2 


P^\s/2 


r^(l+5) 


1/2 


2*4 

f(21) 


Wall  heat  load  per  unit  area, 


^  El _ 1_  2,4 

~  21^  ~  2ttR 


I  /  "a 


'l'  ” 


Wa 


s/2 


(I  +  5) 


1/2 


The  last  two  quantities  are  listed  above  for  completeness,  and 
will  not  be  dealt  with  further  except  to  mention  that  for  rj^  »=  R  the 
maximum  average  enthalpy  per  unit  heat  load  occurs  when  5  =  1;  or 
when  the  power  radiated  is  equal  to  the  power  conducted. 


(22) 
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4.  EXPERIKEOTAL  RESULTS 

A  series  of  measurements  were  made  with  both  argon  and  nitrogen 

in  the  asymptotic  or  fully  developed  discliarge  regime.  The  apparatus 

(4) 

(Fig.  1)  consisted  of  a  "Maecker  type"'  '  stack  of  seven  water  cooled 
copper  segments  of  0.5  cm  thickness,  separated  by  insulating  washers 
of  0.1  cm  thickness.  The  internal  hole  diameter  in  the  stack,  through 
which  the  arc  was  drawn,  was  0.50  cm.  A  special  segment  incorporating  a 
quartz  window  was  constructed  and  installed  midway  in  the  stack  of 
segments  for  measuring  the  power  radiated  per  unit  length  from  the  arc 
column.  Detailed  information  on  the  device  and  the  measurement  tech¬ 
niques  may  be  found  in^^\ 

Arc  characteristic  measurements  (electric  field  as  a  function  of 

current),  taken  at  1.0  atmosphere  pressure,  are  presented  in  Figs,  2 

and  3  with  the  measurements  of  Maecker^ Figures  4  and  5  show  the 

power  radiated  per  unit  length  as  a  function  of  current  at  1.0 

(4) 

atmosphere  pressure,  again  witli  the  measurements  of  Maccker  .  The 
electric  field  as  a  function  of  pressure  for  I  -  200  amperes  current 
is  plotted  on  Figs.  6  and  7  while  Figs.  8  and  9  show  the  power  radiated 
per  unit  length  as  a  function  of  pressure  for  the  same  current. 

A  detailed  comparison  may  now  be  made  between  the  experimental 
measurements  and  the  predictions  of  tfie  simple  theory  (listed  in 
Section  3,  Eqs.  (11)  to  (20)  inclusive). 
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First^  in  order  to  apply  the  theory^  values  for  certain  constants 
must  be  obtained  from  curves  of  arg  .1  and  nitrogen  plasma  properties. 
From  Figs.  10  and  11^  the  following  values  were  determined  for  the 
constants  in  the  linear  theory: 

Argon  Nitrogen 

(fj  =  300  =  6000 


/do\ 

\dl/ 

\  V  a 


1.76 


3.38  X  10' 


fdxx\ 

\d^y 


0.428 


0.720  X  10- 


(23) 


(All  MKS) 

In  addition^  s  was  taken  to  be  zero,  and  t  was  assigned  a  value  of 
1.0  for  lack  of  better  information.  The  predictions  of  the  linear  theory 
may  now  be  calculated  and  drawn  on  Figs.  2  to  9.  The  radius  ratio  at 
first  sight  appears  to  be  implicitly  buried  in  the  equations,  but  may  be 
easily  calculated  and  included  in  the  plots  as  follows: 


1.  Assign  values  to 


s/2 


'da^ 


1/2 


2.  Multiply  each  by  2.4^^^^  from  Eq,  (17),  yielding  El/(l+5) 

3.  Calculate  R/r^^  from  Eq.  (20)  and  5  for  each  pressure  from  Eq.  (11) 

4.  All  quantities  in  Figs.  2  to  9  may  now  be  calculated  from  the 
appropriate  equations.  (For  the  constant  current  plots  of 
Figs.  6,  7,  8,  and  9,  a  simple  cross-plotting  is  needed.) 
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In  general^  it  may  be  seen  that  the  current  and  pressure  dependence 
of  the  radiated  power  from  the  linear  theory  agrees  reasonably  well  with 
experiments  (Figs,  4^  5j  8j  9);  however,  the  prediction  of  the  electric 
field  as  a  function  of  current  does  not  agree  even  qualitatively  with 
the  experiments  as  shown  on  Figs.  2  and  3,  This  disparity  suggests  that 
the  nonlinear  nature  of  the  electrical  conductivity  curve  given  on  Fig.  10 
is  quite  important  for  the  calculation  of  the  electric  field.  Finally, 
after  allowance  is  made  for  the  current  dependence  disparity  in  Figs, 

2  and  3,  the  electric  field  dependence  upon  the  pressure,  presented  in 
Figs,  6  and  7  again  correlates  reasonably  well  between  theory  and 
experiment. 


14 


SIMIIARITY  CONSIDERATIONS 


Up  to  this  pointy  the  linear  theory  has  been  compared  with 

experiments  and  found  to  correlate  reasonably  well  except  for  the 

electric  field  dependence  upon  current.  In  an  attempt  to  account 

for  this  current  dopendenccj  it  is  suggested  that  the  parameter 

in  Eq.  (17)  might  now  be  allowed  to  vary  with  ((|>  -  <{>2^);  because 

from  Eq.  (18),  ((|i  -  (|ip  is  related  to  the  current.  A  similar  process 
_  ,  _  .  .  .  . . ,  . ,  /du\  .  .  . ,  ^  ,  .  ... 


may  be  carried  out  with  the 
forms  are  obtained: 


parameter  and  the  following  similarity 


(1  +  5) 


1/2  R  / 


CO 


Ilii  ( fiV  /  E-V  -  . .  _ ^  Cp  )  /S-V 

|_R(1*5)^'^  (^_)| 

In  order  to  convert  these  expressions  into  forms  in  which 
experimental  results  may  be  conveniently  plotted,  an  expression  is 
needed  for  R/rj^. 

Now,  R/rj^  can  be  written  (see  Eq.  20)  from  the  linear  tlieory  as 


s/2 


(1  +  §) 


R  (1  +  5) 
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Equation  (26)  with  Eq.  (24)  indicates  that 

s/2 


R  (1  +  5) 


1/2 


and  thus  that 


(27) 


and 


(28) 


(29) 


The  similarity  form  of  Eqs.  (28)  and  (29)  may  also  be  obtained  directly 
from  the  energy  equation  using  an  approximation  technique  (see  Appendix 
A  for  details).  All  quantities  can  be  measured  now  and  the  ranges  of 
validity  (if  any)  of  Eqs,  (28)  and  (29)  can  be  checked. 

Using  the  variables  of  Eq,  (28),  Fig,  12  presents  all  of  the 
experimental  data  taken  at  1.0  atmosphere  or  above  in  the  present 
investigation  for  argon,  including  runs  at  different  pressures  and  wall 
radii  (s  is  taken  to  be  zero).  For  some  of  the  operating  conditions, 
the  power  radiated  was  several  times  the  power  conducted  with  large 
increases  in  the  electric  field.  An  empirical  expression  was  used  for 
§  in  order  to  plot  the  R  -  0.377  cm  data  since  radiation  was  not 
measured  at  tl.e  time.  No  data  below  1,0  atmosphere  were  plotted 
because  there  is  some  reason  to  believe  that  Other  effects  become 
important  that  are  not  contained  in  the  simple  theory  upon  which  the 
similarity  parameters  are  based  (non-equilibrium,  for  example),  A 
corresponding  plot  for  nitrogen  is  given  in  Fig.  13. 
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Turning  now  to  Eq.  (29),  Fig.  K  plots  some  experimental  data 
for  argon.  The  pressure  dependence  of  t  1.0  is  satisfactory,  but 
the  absolute  value  of  the  presei  t  measurements  should  not  be  taken  too 
seriously  as  problems  of  window  contamination  were  present  at  tlie  time 
the  data  were  taken. 

On  the  other  hand,  the  measurements  with  nitrogen  on  Fig.  15 
were  taken  after  the  window  contamination  problem  had  been  solved  and 
agree  quite  well  with  those  of  Maecker.  Again  the  pressure  dependence 
of  t  -  1.0  is  satisfactory. 

In  addition  to  using  actual  experimental  results  to  check 
the  similarity  forms  of  Eqs,  (28)  and  (29),  "exact"  calculations 
covering  a  much  wider  range  of  variables  may  be  used.  For  example, 
extensive  calculations  on  nitrogen^^^  are  reduced  using  similarity 
variables  in  Appendix  A,  with  very  convincing  results.  Some  results 
of  this  reduction  are  sketched  in  on  Figs.  13  and  15  for  comparison 
purposes. 

Up  to  this  point  in  the  discussion,  the  question  implicitly 
dealt  with  has  been  "Can  one  collapse  data  taken  over  a  wide  range 
of  operating  conditions  to  a  common  basis?".  The  tentative  results 
of  Figs.  12,  13,  lA,  15,  and  the  plots  of  Appendix  A,  are  quite 
encouraging  in  this  regard  but  await  the  taking  of  considerably  more 
experimental  data. 

From  an  engineering  point  of  view,  one  can  ask  the  following 
question:  "To  what  extent  can  one  predict  the  results  of  operation 
under  a  wide  range  of  conditions  using  simple  formulas,  and  not 
resorting  to  complex  (or  computer)  calculations?".  The  manner  in 
which  the  present  "modified  linear  theory"  may  attempt  an  answer  to 
the  latter  question  is  briefly  discussed  below. 

Assume  that  a  (<|i,  p)  and  u  (^,  p)  are  given.  Essentially, 
in  the  present  formulation,  there  exists  a  correspondence  between 
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the  shapes  of  the  if  curve  in  Eq.  (24)  and  a  P  )j  and  between  the  # 
curve  in  Eq.  (25)  and  u  (ij>,  P^)*  This  correspondence  arises  because 
the  similarity  rules  defined  by  'll'  and  f  were  first  suggested  by  allowing 


and 


the  constant  slopes  of  the  original  linear  approximations^ 

to  now  vary  with  ((^  -  ^j^)  in  the  solutions. 

As  an  example^  if  one  assumes  the  following  forms  for  i|f  and  $: 


*  (x)  =  x" 


(30 


where 


4  (x)  =  x™ 


X  = 


R  (1  +  I) 


1/2 


then  Appendix  B  shows  that 


c  (ij*,  p)  =  Kj^  ((f-  - 


(31 


(32 


and 

u  p)  =  (<f  -  (^^)® 

where  f^  Kj^  g,  and  are  explicitly  given  in  terms  of  n^  m^  and 
C^j  and  vice  versa. 

At  this  point,  we  may  turn  the  procedure  around  and  fit  the 
actual  a  (<f,  p^)  and  u  (^,  p^)  curves  with  the  constants  of  the 
analytical  expressions  for  a  ((|>,  p^)  and  u  ((^,  p^)  (in  the  above 
example,  f,  Kj^,  g,  and  ^2^*  This  process  in  turn  establishes  the 

i(f  and  f  functions,  and  R/t^  by  Eq.  (26). 

Continuing  with  the  example  by  fitting  Eqs.  (32)  and  (33)  to 
argon  and  nitrogen  (see  Figs.  10,  11): 
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Argon 


Ni trogon 


=  500 

=  9000 

f 

=  1/3 

f 

=  1/3 

-  292 

•^1 

210 

g 

=  1/2 

g 

=  2/3 

K2 

3.60  X  10^ 

*^2 

=  2.30  X  10^ 

The  above  values  of  the  fitted  constants  may  be  used  in  turn  in 
the  equations  of  Appendix  B  to  give  values  of  (n,  Cp  m,  C^),  which, 

along  with  Eqs.  (30)  and  (31),  define  the  'I' and  $  functions  for  use  in 
Eqs.  (24)  and  (25).  Equation  (26)  for  R/ r ^  may  then  be  used  to  relate 

the  results  tc  the  experimental  parameters  (Eqs,  (28)  and  (29)  )  and  the 
"modified  linear  theory"  may  be  plotted  on  Figs.  12,  13,  14,  and  15. 

The  modified  linear  theory  may  be  further  reduced,  yielding 
the  ratio  of  power  radiated  to  power  conducted  the  electric  field  E, 
the  arc  current  I,  and  the  power  radiated  per  unit  length 

detailed  in  Appendix  B.  The  results  are  plotted  on  Figs,  2  to  9. 

In  general,  the  proper  trends  with  respect  to  current  and 
pressure  are  exhibited,  and  the  agreement  may  be  said  to  be  "good", 
particularly  considering  the  simplicity  and  resulting  ease  of 
calculation  of  tlie  expressions  used. 
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6.  CONCLUSIONS 

In  order  to  have  a  model  describing  both  the  "heating"  and  the 
"asymptotic"  regions  of  a  wall-stabilized  arc  column^  a  linearized 
analysis  was  made  following  St ine-Watson^  but  including  radiation. 
Attempting  to  account  for  certain  failings  in  the  linear  theory,  the 
solutions  were  "modified"  by  allowing  the  slopes  of  the  linear  ap¬ 
proximations  to  a(ij>)  and  u(<^)  to  vary  witli  (j> ;  the  results  (as  applied 
to  the  asymptotic  region.  Initially)  are  as  follows: 

1.  The  functional  dependences  of  the  various  parameters  are 
now  approximately  correct  compared  to  the  linear  theory 
(particularly  the  electric  fie  d  as  a  function  of  current). 

2.  The  curve  fitting  of  the  plasma  property  curves  ^(il*)  and 
u(^)  usi  ng  simple  power  laws  is  considerably  more  flexible 
and  less  arbitrary  than  ti>c  straight  line  approximations 
in  the  linear  theory.  The  fact  that  o(<j>)  and  u(<|>)  shared 
a  common  constant  parameter  did  not  seem  to  be  too 
restrictive. 

3.  An  "effective"  arc-to-wall  radius  ratio,  ^  ^  is  defined 
wJiich  may  be  followed  through  the  equations  and  studied  in 
the  solutions. 

A,  Similarity  forms  are  suggested  by  the  modified  linear  theory 
or  alternately  tlirough  approximations  to  the  basic  equations 
these  are  quite  useful  for  plotting  exper imenua 1  results  and 
collapsing  theoretical  calculations,  and  afford  a  direct 
comparison  between  the  two. 

5.  For  scaling  and  predicting  purposes,  the  results  are 

presented  in  alternate  forms  such  that  the  best  information, 
whether  experimental  o-  theoretical,  may  be  used.  (For  exam 


experimental  curves  at  1.0  atmosphere  may  be  scaled  as  a 
function  of  pressure  using  theoretical  powers  of  pressure, 
and  analytical  expressions), 

6.  There  is  strong  evidence  tliat  tlie  "deviations  from  experiment 
of  an  exact  solution  using  present  "state  of  the  art"  calcu¬ 
lations  of  plasma  properties  are  at  least  as  large  as  the 
deviations  introduced  by  the  approximations  in  the  modified 
linear  theory.  Tlie  fact  that  functions  are  known  for  argon 
and  nitrogen  plasma  properties  that  will  reproduce  the 
experimental  results  (at  atmosplieric  pressures)  in  an  exact 
calculation  is  beside  the  point  because  they  were  derived 
precisely  with  that  constraint.  Hence,  in  a  practical  sense, 
there  is  little  value  in  using  exact  numerical  calculations 
for  arc  column  analysis,  even  for  the  prediction  of  absolute 
magnitudes.  An  interesting  trend  noted  is  that  the  modified 
linear  theory  and  tlie  AVCO^^^  plasma  properties  seem  to  have 
their  "errors"  in  opposite  directions  so  that  using  approxi¬ 
mations  to  the  AVCO  plasma  properties  in  the  mc?dified  linear 
theory  gives  better  agreement  than  tlie  "exact"  calculations. 
The  above  chance  effect  was  noted  with  nitrogen  and  with 
argon  (using  AVCO-type  argon  calculations)  but  might  be 
expected  to  hold  with  other  gases,  since  similar  fcniiulas 
and  shapes  of  curves  are  involved. 

Eventually  the  real  value  of  such  methods  as  the  modified  linear  theory 
may  be  shown  in  the  gas  heating  region  by  providing  the  insight  needed 
to  incorporate  additional  items  (such  as  radiation  and  arc  radius)  into 
the  simple  tlicorics. 

Experiments  in  the  gas  heating  region,  including  the  measurement 
of  radiation,  are  urgently  needed  for  further  comparisons. 
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APPENDIX  A 

DERIVATION  OF  THE  ASYMPTOTIC  REGION  SIMILARITY 
RELATIONSHIPS  FROM  BASIC  EQUATIONS 


The  well-known  Elenbaas-He Her  differential  equation  for  an 
optically  thin^  steady  state^  axially  synunetric  arc  of  radius  R  and 
constant  electric  field  strength  E  may  be  written  in  the  following 
form  (using  Ohm’s  Law): 


u  P) 
c  (<l>,  P)  E^-J 


( 


at 


0 


at  p  ««  1^  o  *»  0 


In  the  above  equation. 


fi  «  r/R 


i 

^  ^  K  (T,  p)  dl  =  ^  (Gas,  I,  p,  R,  p,  T^) 


w 


E  =  E  (Gas,  I,  p,  R,  T^) 


and  T^  is  the  wall  temperature  (r-  0). 


Assui.je 


i)  The  ratio  of  u  to  a  is  a  function  of  pressure  only. 


P) 

o  (<?/  P) 


F  (P) 


(3 
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and 


ii)  The  ratio  of  the  volumetric  radiated  power  to  the 

volumetric  input  power  is  the  same  as  the  ratio  of  the 
power  radiated  per  unit  length  to  the  power  input  per 
unit  length. 


P) 

a  (4^  p) 


=  (i^y  i) 

E" 


Under  the  assumption  inherent  in  Eq.  (34)^ 


(36) 


El  -  P  .  +  P  . 

rad  cond 


(37) 


where  P  ,  is  the  power  conducted  to  the  wall  per  unit  leneth. 

cond  I  • 

Then 


1  - 


--H-lii  p) 

O  (<j>,  p)  E^J 


rad 


P  j  +  P 
rad  cond 


1 


where 


1  +  I 
P 


I  -  ?  (Gas,  I,  p,  R,  T^)  .  ^ 

rad  cond 


Eq,  (34)  becomes 


L(1  +  5) 

The  total  current  may  be  obtained  from  Olim's  I.aw 


R 


I=^j2'^rdr  =  R  E^c(<j' 


,  p)  2-  Pd  p 


(38) 


(39) 


(40) 


Dividing  both  sides  by  R  (1  and  considerirtg 


(f,  p)  -  c 


ER 


.(1  -  ?) 

the  followinc  form  results; 


1/2 


,  P,  C  from  Eq.  (39), 


I 


R  (1 


fn 


ER 


(1  ^ 


9  J  e 


(4n 
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Turning  to  the  radiation  now,  from  Eq.  (36) 


^.d  ‘  E 


Using  Eq.  (Al), 


(1  +  5) 


R  F  (p) 


(1  +  5) 


The  pressure  dependence  in  Eqs,  (Al)  and  (43)  can  be  completely 
included  in  the  similarity  expressions  if  the  following  further  assump- 


t ions  are  made : 


o  (<[•,  P  )  =  a  (<^,  p^)  I  ^ 


and  iv) 


wliere  p  is  a  reference  pressure  (taken  to  be  atmospheric  for  convenience), 

3 

/u\  /u\ 

and  (— >  =  — TT - r  is  some  average  value  of  /  —  i  over  the  range  of 

V/a  <7  (f,  p^)  la/ 

interest,  evaluated  at  the  reference  pressure. 

These  further  assumptions  may  be  carried  through  the  equations  and 

tlie  final  results  summarized  as  follows: 


o  (f,  P)  -  c  (<j>,  p  ) 


a'  p 


„  © 


a  (f,  p)  e' 


a  /£_Y  ^  _  ^rad 

e'  Ua  ^ 
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then 


ER 


8/2 


(1  +  5) 


1/2 


R(1  +  ?) 


1/2 


(48) 


(49) 


(50) 


(51) 


rp  ,  E 

rad 


(52) 


I  - 


Of  the  quantities  in  the  above  equations^  1,  E,  ^  cad 

P  ,  .  .  ^ 

are  all  measurable  and  have  definite  val  ’.es  at  each 


El  -  P 


exper imental^operat ing  point. 

The  quantities  s  and  t  are  theoretical^  or,  alternately,  fitted 
from  experiments. 

The  quantity  may  be  obtained  either  theoretically  from  Eq. 

(51),  or  experimentally  from  Fq.  (52)  by  averaging  experimental  quan¬ 
tities  over  the  region  of  interest  at  the  reference  pressure  p^. 

An  interesting  check  of  the  above  similarity  forms  is  provided  by 
plotting  the  results  of  extensive  arc  column  calculations  using  the 
above  similarity  variables;  for  example,  the  calculated  "exact"  nitrogen 
solutions  of  which  cover  a  very  wide  range  of  parameters,  may  be 
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so  plotted. 

First,  the  usual  non-radiating  similarity  parameters  ER  vs  l/R  are 

plotted  In  Fig.  16.  Each  individual  symbol  represents  one  calculated 

cond^'tion  obtained  by  assuming  values  for  the  electric  field  and  axis 

temperature  and  then  calculating  the  other  parameters  of  interest 

(including  wall  radius).  In  general,  calculations  of  this  nature  are 

difficult  to  compare  directly  with  experiment  because  of  the  varying 

wall  radius;  however,  if  the  similarity  expressions  hold,  this  variation 

(4) 

will  be  of  no  consequence.  The  measurements  of  Maecker  are  also 
dashed-in  on  Fig.  16  for  reference.  The  effect  of  radiation  in  this 
plot  is  to  scatter  points  upward  from  base  lines  of  zero  radiation  at 
each  pressure. 

Next,  the  similarity  parameters  for  constant  pressure  are  plotted 
in  Fig.  17.  The  use  of  these  parameters  very  effectively  collapses 
most  of  the  points  to  the  appropriate  zero-radiation  curve  for  each 


pressure. 

Now  an  attempt  will  be  made  to  further  reduce  the  calculations  by 
including  pressure  explicitly  in  the  plotting  parameters.  Ref.  1  gives 
fitted  expressions  for  a,  (Ji,  and  u  which  may  be  applied  to  nitrogen 
above  30,000  °K  (See  Zqs.  (65),  (67)  and  (73)  of  Ref.  1).  From  these 

equations,  -  ^  is  insensitive  to  .Ji  (fulfilling  assumption  i). 


equations,  —  ^  is  insensitive  to  .ji  (fulfilling  assumption  i), 

s  =  0.154  and  t  =  1.806.  Accordingly,  Figs.  18  and  19  may  be  plotted. 

As  expected,  the  pressure  dependence  over-compensated  at  the  low  temper¬ 
ature  end  of  Fig.  18  but  in  general  did  a  remarkable  job  of  collapsing 
the  data  in  the  two  figures.  Particularly  satisfying  in  Fig.  19  was  the 

l*rad  I 

manner  in  which  (compared  to  a  — —  vs  —  plot)  tht  5  parameter  adjusted 

R^  ^ 

the  points  at  each  pressure  to  fall  on  continuous  similar  curves,  the 

/Pa\ 

Parameter  moved  the  (similar)  curves  at  the  higher  pressures  to 


the  left,  and  the 


c-y 


parameter  collapsed  the  curves  downwards  to 


the  atmospheric  pressure  cu'.ve.  Figs.  17,  18  and  19  also  indicate 
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that  s  and  t  values  can  probably  be  assigned  such  thoL  the 
similarity  forms  may  be  used,  at  least  piecewise,  over  the  complete 
region  of  interest. 

Finally,  a  direct  comparison  may  now  be  made  between  the  calcu- 
lations^^^  and  the  experiments  of  Maecker^^^  for  nitrogen  with  no 
problems  involving  different  operating  conditions  confusing  the  results. 
The  disagreements  shown  on  Figs.  17,  18,  19,  or  Figs.  13,  15  occur 
essentially  because  the  AVCO^^^  electrical  conductivity  is  as  much  as 
607,,  higher  and  the  volumetric  radiated  power  as  much  as  a  factor  of  5 
higher  than  the  equivalent  properties  obtained  from  the  inversion  of 
Maecker's  measurements  in^^\ 
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APPENDIX  B 


THE  MODIFIED  LINEAR  THEORY  SOLUTION 


From  the  exact  linear  solution,  Eqs.  (17),  (18)/  and  (12)  combined 
with  Eq.  (11),  may  be  written  as  follows: 


(53) 


(54) 


(55) 


and 


vd^4^ 


the  "average"  slopes  of  a(<|>,  p  )  and  u(i|>,  p  )  respectively,  to  become 


functions  of 


)  in  the  solution. 


Considering  and  to  be  functions  of  (f-'j'i^)/  then  from 


</  a 


Eqs.  (53),  (54),  and  (55) 


R  (1  ^ 


(56) 
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and 


To  simplify^  assume  the  following  forms  for  ij(  and  f: 


♦  (X)  -  x" 

(58) 

i  (x)  =  C2  x"* 

(59) 

where 


Then^  from  Eqs.  (53)^  (56),  and  (58), 


-n 

X 


(from  Eq.  54) 


(60) 


Also,  from  Eqs.  (55),  (57),  and  (59), 
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Substituting 


F  ina  1  ly 


where 


a  may  be  evaluated  in  terms  of 

2 


Pa>  ^  ^  ^ - =rj^ 

™«  [2-TJ,  (2.4)3^ 

(62) 

-  Kj  (.f  - 

(63) 

£  .  i2Jl 

l+n 

(64) 

0 

(z.A^ 

(65) 
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or 


n  = 


1-f 

1+f 


C,  = 


2.4 


1  Ip  — iliJL 

Ki  “  12^  Ji(2.4)j  1+f 


Similarly^  u  may  be  evaluated  in  terms  of 

^m-n- 1  ^ 

I —  “n\  n  + 1  / 

I^TT  (f  .fp 


Pa> 


m 

,  _  2  4  C 

.  P  ■  Vl  I  L 


m 

n-1 


E  K2 


where 


m 

n+1 


=  2.4  C 


m 

n+1 


2  I  C 


2Tr 


•Ij  (2.4)1 


fer) 


or 


m 


_  2£_ 
1+f 


C.  = 


2  -  2.4  ~zr^ 

1+f  2n  (2.4)  '  1+f  / 


(66) 


(67) 


(68) 

(69) 

(70) 

(71) 

(72) 

(73) 
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Hence^  If  che  group  of  parameters  (n^  m,  C2)  is  known^  the  group 
(f,  g,  K^)  may  be  calculated  from  Eqs.  (64),  (65),  (/O),  and  (71). 
Conversely,  if  the  group  (f,  Kj^,  g,  K2)  is  known  (for  example,  by 
curve  fitting  0(1^,  p^)  and  u(^,  P^),  the  group  (n,  Cj^,  m,  C^)  may  be 
calculated  from  Eqs.  (66),  (67),  (72),  and  (73). 

For  example,  assume  that  the  experimental  results  can  be  plotted 
using  the  following  similarity  forms  (including  fitting  values  of  s 
and  t,  if  needed) ; 


If  a  reasonable  value  of  can  be  assigned  (say  from  fitting  theoretical 
curves  of  o(<{>,  p^)  and  u(|>,  p^)  near  the  axis),  then  from  Eqs.  (26)  and 
(74) 


(76) 


and  the  similarity  forms  in  Eqs,  (56)  and  (57)  may  be  plotted.  These 
plots  may  bo  fitted  witii  (n,  Cj^,  m,  C^)  from  Eqs.  (58)  and  (59)  and 
thus  (f,  Kj^,  g,  K^)  may  be  calculated.  The  shapes  of  the  a(i}>,  p)  and 
u(j',  p)  curves  implied  by  the  experimental  results  arc  then  given  by: 
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s 


(77) 


a((Ji,  p)  = 


(78) 

Conversely,  the  similarity  functions  of  Eqs.  (74)  and  (75)  may  be 

obtained  from  the  fitted  shapes  of  the  a(<i>,  p  )  and  u(i,  p  )  curves 

(yielding  f,  g,  K^)*  After  (n,  Cj^,  m,  C^)  is  calculated,  then 

~  may  be  calculated  as  a  function  of  x  from  the  following  equation 
’^l 

(see  Eqs.  (26),  (56),  and  (58)  ): 


^1 


exp 


X 


n+1 


(79) 


where 


R  R 

Knowing  —  =  —  (x),  then  the  similarity  functions  of  Eqs.  (74)  and 
*^1  ^1 

(75)  immediately  follow  using  Eqs.  (58)  and  (59). 

If  it  is  desired  to  obtain  such  quantities  as  electric  field, 
current  and  radiated  power  per  unit  lengtli  for  a  specific  geometry, 
an  expression  is  needed  for  5  from  the  modified  theory. 

The  original  expression  for  |  was  the  following:  (See  Eq.  11) 
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from  equations  of  the  modi  fil’d  linear  theo-y,  with 


hence  after  values  of  s  and  t  are  detennined  in  addition  to 
f,  Kj,  g  and  K^)  by  curve  fitting  o(i|),  p)  and  u(ij),  p)  with  the  ex¬ 
pressions  of  Eqs.  (77)  and  (78),  all  quantities  of  interest  may  be 
obtained  as  a  function  of  x  for  givim  values  of  R  and 


If  constant  current  calculations  are  desired,  simple  cross 
plotting  is  necessary  (as  with  tlu*  lini’ar  theo  ry). 

Accordingly,  values  of  §  have  been  calculated  as  functions  of 
current  and  pressure  for  argon  and  nitrogi’n.  The  results  arc  plotted 
in  Figs.  20  to  23  along  with  experimental  values  and  the  predictions 
of  the  straight  linear  theory.  The  re  la  t  i  onsii  i  ps  between  the  irodified 
linear  theory  and  the  exp'Tinents  ari’  wtiat  one  would  expect  after 
examining  the  curve  fits  i>n  Fig.  11.  Figs.  20  to  23  may  be  used  in 
turn  to  reduce  the  similarity  form  of  the  modifii’d  linear  thr-ory  given 
on  Figs.  12  to  15  back  to  th.  quantities  plotted  in  Figs.  2  to  9  for  a 
specific  geometry.  The  res'ults  are  discussed  in  Section  5. 

The  natural  question  arisi’S  as  t’O  thf-  relation  between  the 
expression  for  5  given  by  Eq.  (80)  and  rl;at  given  by  Eq.  (50)  in 


(8( 
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Appendix  A,  which  was  derived  using  assumptions  in  the  basic  equations. 
Applying  the  assumptions  of  Appendix  A  to  the  present  theory, 

f  "  8  (81) 


and 


From  Eqs.  (64)  and  (70), 


Hence 


or 


(S?) 


(83) 


(84) 


From  Eqs.  (67)  and  (73), 


(from  Eq.  82). 


ERI  ^ 


s/; 


After  replacing  — C  x  with  - 

"l  (1*5)*^^ 


from  the  present  t'neory, 


Eq.  (84)  becomes  identical  with  Eq.  (^0)  combined  with  Eq.  (48)  of 
Appendix  A. 
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Honcc  the  modified  linear  theory  reduces  to  the  results  given  in 
Appendix  A  when  o(ij»,  p)  and  u(i^,  p)  are  required  to  have  the  same  ^ 
dependence . 
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DETAIL  OF  ARC 
COLUMN 

NOTE: 

1.  ALL  COMPONENTS  WATER  COOLED  UNLESS  OTHERWISE  NOTED 

2.  ALL  INSULATORS  BETWEEN  COMPONENTS  ARE  A  COMPOSITE  OF 
BORON  NITRIDE  ,  TRANSITE  8  SILICON  RUBBER 

FIC.  1  ARC  COLUMiN  SCHKM^XTIC 
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ARGON  GAS 
PsI.O  ATM 


(J3iauj/S|I0A  ^0I)‘3  an3ld  Dldl03“l3 
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FIG.  2  ELECTRIC  FIELD  AS  A  FUNCTION  OF  CURRENT  -  ARGON  GAS 


MODIFIED  LINEAR  THEORY 


FIG.  3  ELECTRIC  FIELD  AS  A  FUNCTION  OF  CURRENT  -  NITROGEN  GAS 


ARGON  GAS 
P=I.O  ATM 
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FIG.  4  POWER  RADIATED  PER  UNIT  LENGTH  AS  A  FUNCTION  OF  CURRENT  -  ARGON  GAS 
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Flu.  5  POWER  RADIATED  PER  UNIT  LENGTH  AS  A  FUNCTION  OF  CURRENT  -  NITROGEN  GAS 


ARGON  GAS 
1=  200  AMPS 
R  =  0.25cfn 
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ARC  PRESSURE  P/Pq  (atmospheres) 

FIG.  6  ELECTRIC  FIELD  AS  A  FL’NCTION  OF  PRESSURE  -  ARGON  GAS 


MODIFIED  LINEAR  THITORY 
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THERMAL  CONDUCTIVITY  INTEGRAL  ( 10^  watts/meter) 

FIC.  10  ELECTRICAL  CuNDUCTlVlTY  AS  A  FUNCTION  OF  IP  THER>L'L  CONUUCTIVI TY  INTEGRAL  FOR  ARGON 


U  =  338*lO  (6-300) 
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FIG.  IJ  arc  CH.'J  ACTKRI.STJC  V.AI.GKS  r-OTI'En  USING  SIMII^RITY  VARIABIFS 
NlTRlV.KN  GAS 
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FIG.  14  RADIATED  POWER  PER  UNIT  LENGTH  PLOTTED  USING  SIMILARITY  VARIABLES 
ARGON  GAS 
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FIG.  16  C.MXULATED  NITROGEN  CVL\RACTER1STICS  PLOTTED  WITH  (N  %  .VITNG)  SIMILARITY 

variables 
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INCLUDING  RADIATION  AND  PRESSURE 


EXPERIMENTS)  ARGON  GAS 


ooooooooo 
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}  diMOd  aBionoNOO  oi  asiviavd  do  oiiva 
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FIG.  21  RATIO  OF  RADIATED  TO  CONDUCTED  POWER  AS  A  FUNCTION  OF  CLTRRENT  -  NITROC^EN  GAS 
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FIG.  22  RA'.’IO  OF  RADIATED  TO  CONDUCTED  POWER  AS  A  FUNCllON  OF  PRESSURE  -  ARGON  GAS 
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FIG.  23  RATIO  OF  RADIATED  TO  CONDUCTED  POWER  AS  A  FUNCTION  OF  PRESSURE  -  NITROGEN  GAS 


